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Abstract: Calotropis procera (Aiton) Dryand (Apocynaceae) is a native species in tropical and subtropical 
Africa and Asia. However, due to its fast growing and drought-tolerant, it has become an invasive species 
when it was introduced into Central and South America, as well as the Caribbean Islands. Currently, C. 
procera displays a wide distribution in the world. Invasiveness is important, in particular, because many 
invasive species exert a high reproductive pressure on the invaded communities or are highly productive in 
their new distributed areas. It has been suggested that a very deep root system and a high capacity to 
reduce stomatal conductance during water shortage could allow this species to maintain the water status 
required for a normal function. However, the true mechanism behind the successful distribution of C. 
procera across wet and dry environments is still unknown. C. procera leaves were collected from 12 natural 
populations in Brazil, Colombia and Mexico, ranging from wet to dry environments during 2014—2015. 
Many traits of morphology and anatomy from these distinct morphotypes were evaluated. We found that 
C. procera leaves had a considerable capacity to adjust their morphological, anatomical and physiological 
traits to different environments. The magnitude of acclimation responses, i.e., plasticity, had been 
hypothesized to reflect the specialized adaptation of plant species to a particular environment. However, 
allometric models for leaf area (LA) estimation cannot be grouped as a single model. Leaves are narrower 
and thicker with low amounts of air spaces inside the leaf parenchyma in wet environments, while they are 
broader and thinner with a small number of palisade cell layers in dry environments. Based on these, we 
argue that broader and thinner leaves of C. procera dissipate incident energy at the expense of a higher rate 
of transpiration to survive in environments in which water is the most limiting factor and to compete in 
favorable wet environments. 
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1 Introduction 


Invasive plants are non-native species that affect the biodiversity of the regions into which they 
are introduced or that they colonize (Colautti et al., 2004). In the last 20 years, local invasions by 
non-native species have become a focus of ecological research, primarily because of the extensive 
impact they have on ecosystems or the set of new ecological strategies they can play in colonizing 
new environments at a small geographical scale or at the macroscale (Alberio and Comparatore, 
2014; Pagnucco et al., 2015; Dawson et al., 2017). Invasiveness is important, in particular, 
because many invasive species exert a high reproductive pressure on the invaded communities 
(Simberloff et al., 2012) or are highly productive in their new distributed areas (Blossey and 
Notzold, 1995). Calotropis procera (Aiton) Dryand (Apocynaceae), known as milkweed, is an 
example of an invasive species (Muriira et al., 2018). It is native to tropical and subtropical Africa 
and Asia and is introduced into Central and South America, as well as the Caribbean Islands. 
Currently, C. procera displays a wide distribution (Singh, 1995) ranging from dry environments 
to wet environments, even grows in the area where the elevation is 2000 m a.s.l. (Hassan et al., 
2015). C. procera can quickly establish itself as a weed along degraded roadsides, lagoon edges 
and overgrazed native pastures. Morphological variability exists among different provenances 
(populations) of this species (Peguero-Pina et al., 2014). 

C. procera is an evergreen and shrubby plant with the ability to maintain high rates of net 
photosynthesis, even during the dry season (Rivas et al., 2017). This ability is a key feature that 
contributes to the competitive success of C. procera under diverse environmental conditions. 
Ecophysiological adaptations, such as the presence of a very deep and stout taproot reaching 
depths ranging from 1.7 to 3.0 m in Indian sandy desert soils may aid in plant survival in extreme 
drought environments (Sharma, 1968). In addition, low stomatal conductance during water 
shortage has been ascribed as an important mechanism to maintain relatively high water content 
and water potential that described these phenomena as the intrinsic water use efficiency (WUE) 
(Rivas et al., 2017), forming this species a classical desiccation postponement morphotype 
(Antunes et al., 2017). 

Qualifying plant form and function at a large scale is the first step in improving the predictive 
capability of vegetation models (Wang et al., 2013; Diaz et al., 2016). It is well known that leaves 
are highly sensitive to environments (Fahn, 1990; Wright et al., 2004; Evert, 2013). The internal 
structure of leaves, such as the arrangement of palisade and spongy parenchyma, affects the 
diffusion of gases and partially determines the photosynthetic and non-photosynthetic tissues 
(Peguero-Pina et al., 2016a). The specific leaf area (SLA, ratio of leaf area to leaf dry weight) is an 
indicator of the morphological and anatomical traits (Villar et al., 2013). In general, low SLA is 
mainly due to the thicker leaf caused by a large volume of mesophyll and air spaces, as well as a 
high leaf density due to the high volumetric fractions of mesophyll and low fractions of air spaces 
and epidermis (Villar et al., 2013). Decreases in SLA have been postulated as the results of denser 
leaf laminae, associated with lignification and thickness of cell walls cuticle waxes, trichomes and 
other non-photosynthetic structures, here referred to as sclerophylly (Villar et al., 2013). On the 
other hand, low SLA is often considered as an indicator for high photosynthetic potential (Reich 
et al., 1992). Poorter et al. (2006) found that low SLA is related to a long leaf life span and a 
conservative growth rate, with higher construction costs for the leaf area. In addition, thicker 
leaves display a greater light absorbance (Falcioni et al., 2017), in part due to chloroplasts stacked 
in a small area, while thinner leaves spread chloroplasts throughout a larger leaf area. Both 
physiological strategies are common in response to light intensity (Evans and Poorter, 2001; 
Poorter et al., 2009). Large leaves have smaller convection coefficients and a higher resistance to 
heat transfer through leaf boundary layers than smaller leaves, and leaf size may change in order 
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to optimize leaf temperature (Leigh et al., 2016). Larger leaves may be intrinsically vulnerable to 
drought-induced embolism due to the lower vein length and larger xylem conduit diameters 
(Peguero-Pina et al., 2014; Gil-Pelegrin et al., 2017). These findings could be extrapolated to the 
leaves that develop in arid (like Mediterranean) environments (Peguero-Pina et al., 2014; 
Peguero-Pina et al., 2016b). Sclerophylly is perhaps an important mechanism to adapt to the 
extremely drought and hot environments with high light intensities. However, the mechanisms 
facilitating successful distribution of C. procera across wet and arid environments are still 
unknown. 

The aim of the study was to determine the morphological traits of C. procera leaves that allow 
it successfully distributed in wet or dry environments and become an invasive species. In addition, 
we evaluated which anatomical and morphological trait of leaves is related to drought. For this 
purpose, we hypothesize that C. procera had various anatomical and morphological traits 
according to their provenances. 


2 Materials and methods 


2.1 Plant materials 


C. procera leaves were collected from 12 natural populations along a macro-environmental 
gradient (i.e., rainforest, arid and semi-arid regions) in Mexico, Colombia and Brazil during 
2014-2015. About 10 populations were from Brazil, 1 from the Caribbean Atlantic Region 
located in Colombia, and 1 from rainforest located in the mountain of southern Mexico (Fig. 1). 
Geoclimatic details are described in Table 1. Climatic datasets from 12 populations during last 15 
years were downloaded from https://pt.climate-data.org. The aridity index was estimated by 
Martonne aridity index (MAI; MAI=P/(T+10), where P is the annual precipitation (mm) and T is 
the monthly mean temperature (°C)). Based on the MAI and climatic characteristics 
(Sant'anna-Neto et al., 2015), we classified Brazilian populations as: (i) Atlantic Rainforest 
(MAI>35), (ii) Semi-arid region (25<MAI§<35), and (iii) Arid region (MAI<25). The population 
collected in Montería, Colombia is referred to as Caribbean and the population collected in 
Oaxaca, Mexico is referred to as Mountain Rainforest (Fig. 1). 
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Fig. 1 Population distributions of Calotropis procera in Brazil, Colombia and Mexico. The states of Alagoas, 
Pernambuco, Paraiba and Minas Gerais in Brazil are shown. The mean annual precipitation during the last 15 
years in all populations is displayed. Different lowercase letters indicate significance of the mean annual 
precipitation among different populations at P<0.05 level. 


2.2 Leaf morphology and specific leaf area (SLA) measurements 


A total of 5937 fully expanded leaves were harvested (approximately 500 leaves of each 
population) in the field at the end of growing season during 2014—2015. The leaves were then 
scanned at a 1200x1200 dpi resolution (HP, PSC1410, Palo Alto, CA, USA) and the images were 
analyzed using the Image Pro® Plus software (version 4.5.0.29, Media Cybernetics, Silver Spring, 
MD, USA). The maximum leaf length (L; from the lamina tip to the intersection point of petiole 
and midrib), the maximum leaf width (W; the longest distance perpendicular to midrib) and 
observed leaf area were measured by digitizing the leaf contours with the aforementioned image 
software. To measure SLA (m7/kg), we oven-dried 50 randomly selected leaf discs (12 mm in 
diameter) from each population for 72 h at 70°C and weighed (Witkowski and Lamont, 1991). 


2.3  Allometric model for leaf area estimation 


Allometric models were run following the methods described in Antunes et al. (2008) and 
Pompelli et al. (2012) for estimating leaf area (LA) using L, W or LxW. The estimated leaf area 
was taken as dependent variables and observed leaf area as an independent variable (n=500). 
After an initial run, the preliminary models were evaluated for bias. The best unbiased models 
were validated using an independent data set (n=200) from the same morphotype. The validated 
and unbiased estimated equation (power models) from each morphotype was then tested for 
equality by applying the procedure described in Graybill (2000). In sum, the first step of this 
procedure was tested if all equations are equal (a<0.05; F test). If not, the population was 
excluded and another was pooled in a single model. This procedure was repeated for all possible 
combinations until a>0.05 (F test). 


2.4 Leaf anatomy and leaf symmetry measurements 


Leaf anatomy was compared among the following four populations that were highly divergent in 
terms of leaf morphology (i.e., SLA and L/W ratio), mean annual precipitation and MAI: (i) in the 
Atlantic Rainforest of the Pernambuco State (PE) (Köppen climatic classfication: equatorial 
climate)—Paulista-PE (#3), (ii) in the arid region of the Pernambuco State (Köppen climatic 
classfication: desertic climate)—Afranio-PE (#10), (iii) in the Caribbean region, where the climate 
is strongly influenced by North Pacific ocean currents (Köppen climatic classfication: 
semi-humid tropical climate)—Monteria-CO (Colombia) (#11), and (iv) in a mountain rainforest 
with annual precipitation of about 2170 mm in Mexico (MX) (Köppen climatic classfication: 
semi-humid tropical climate)-Oaxaca-MX (#12) (Table 1). All sampled tissue were fixed 


Table 1 Locations and basic characteristics of Calotropis procera from 12 populations 


Population Ecosystem/ Elevation Mean annual 


node City-state/country Country Latitude Longitude provenance (m) precipitation (mm)! MAI 
#01 Aimorés-MG Brazil 19°29'33"S 41°03'59"W Atlantic 80 1115+81° 37.2 
#02 Jaboatão Guararapes-PE Brazil 08°11'60"S 34°55'60"WŴW Atlantic 6 1811493" 53.0 
#03 Paulista-PE Brazil 07°51'02"S_ — 34°50'54"W Atlantic 4 1883493" 45.6 
#04 Pombos-PE Brazil 08°08'48"S = 35°23'35"W Semi-arid 237 1066+56° 31.2 
#05 Santa Terezinha-PB Brazil 06°59'13"S 37°18'08"W Semi-arid 270 758468% 28.6 
#06 Alagoinha-PE Brazil 08°28'40"S 36°47'23"WŴW Semi-arid 690 636+66% 25.4 
#07 Jaramataia-AL Brazil 09°39'01"S 37°00'22"WŴW Semi-arid 176 7434444 34.9 
#08 Serra Talhada-PE Brazil 07°54'35"S_ — 38°17'59"W_ Arid region 429 651451% 24.8 
#09 Parnamirim-PE Brazil 08°04'21"S 39°34'00"WŴ Arid region 412 478+46° 15.6 
#10 Afrânio-PE Brazil 08°29'12"S 41°03'21"WŴ Arid region 579 474+49° 13.1 
#11 Monteria-CO Colombia 08°25'47"N_ — 75°53'24"W Caribbean 18 1380+65° 40.6 
#12 Oaxaca-MX Mexico 17°04'45"N 96°46'24"W Eta 2360 2166+83* 72.2 


Note: ! Climatic data in all 12 populations during the last 15 years were downloaded from https://pt.climate-data.org. Different lowercase letters 
indicate significant differences among different populations at P<0.05 level. Mean+SE; n=15. MAI, Martonne aridity index; MG, Minas Gerais; 
PE, Pernambuco State; PB, Paraiba; AL, Alagoas; CO, Colombia; MX, Mexico. 
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formalin-acetic-alcohol (FAA; 10% formalde-hyde, 5% acetic acid and 50% ethanol), dehydrated 
(ethanol series) and embedded in methacrylate. Paradermic sections (4 um) and cross-sections (7 
um) were stained with toluidine blue. Anatomical data were quantified using the Image Pro® Plus 
software. The following anatomical data were measured: (i) total leaf thickness, (ii) adaxial and 
abaxial epidermal thicknesses, (iii) palisade and spongy parenchyma thicknesses, (iv) adaxial 
palisade number, and (v) palisade and spongy air spaces, computed as a percentage of total 
cross-sectional area of the mesophyll tissue. With the exception of palisade and spongy air space, 
all measurements were conducted randomly 10 times per micrograph, using 20 light micrographs 
per morphotype. To better illustrate the differences between morphotypes, we traced drawings 
from photomicrographs using a representative sample. 

For leaf symmetry measurements, 50 fully expanded leaves per morphotype were measured. 
Two different morphological traits of the leaves were directly measured: (i) leaf width and (ii) 
absolute interrib distance, using the method described by M@ller and Eriksson (1994). Data of 
measured leaf characters were tested by Shapiro-Wilk test. When a non-normal distribution was 
observed, the values were normalized by applying the Box-Cox transformation, as recommended 
by Souza et al. (2005). 


2.5 Plasticity index 


Plasticity index was associated with the analysed morphological and anatomical traits, ranging 
from 0.0 (no plasticity) to 1.0 (maximal plasticity). The index was calculated as the difference 
between the minimum and the maximum values divided by the maximum average value of each 
leaf morphological and anatomical trait (Valladares et al., 2002). 


2.6 Pearson's correlation coefficient and principal component analysis (PCA) 


The Pearson's correlation coefficient was used to examine the relationships among variables. 
Similarity between groups of data, were calculated using an adjusted r coefficient through 
Corrplot Package in R software (Wei and Simko, 2016). To test whether the alteration patterns in 
specific sets of traits (morphological and anatomical) were consistent among different populations, 
we used data from all 12 populations of C. procera. Magnitude problems were avoided by 
standardize all variables before the main PCA components (Z-score), with the goals of controlling 
the differences in each measurement variable and also comparing measurements that did not 
present the same units (Gotelli and Ellison, 2012). For this, we used Z standardization, thereby 
satisfying all the basic requirements for testing hypotheses, including data normality and 
homogeneity of variance (Gotelli and Ellison, 2012). Different analyses of variance (one-way) 
were complemented by Tukey's test, and the full range of traits (i.e., multivariate attribute vector) 
associated with different populations were compared. For this, we used the most variable axes 
scores of the data (i.e., PC1 (principal component 1) and PC2 (principal component 2)) (Gotelli 
and Ellison, 2012). This procedure was performed for both morphological and anatomical traits. 
Finally, we gathered all the layers of the PCA on a standard scale, joining them in a single graph 
with the aid of the ggplot2 and Vegan packages in R software (Husson et al., 2017). 


3 Results 


3.1 Allometric model for leaf area estimation 


Allometric models were used to predict leaf area from their linear dimensions (L, W or LxW). The 
best unbiased equation was a power model: Y=,(L x W)"e,, where Bo and Bı are the coefficients; 
L is the leaf length; W is the leaf width; and e; is the random error that represents the discrepancy 
in the approximation (Fig. S1). The equation was validated and the estimated parameters (B's) for 
model calibration were not statistically different (P>0.05) from the validating equations using 
independent sets of leaves from the same morphotype. The majority of the calibrating equations 
estimated for a single morphotype were not statistically similar (P>0.05) to each other by the 
model identity test (Table 2). Thus, we grouped the morphotypes as the order of G1: #12; G2: #10; 
G3: #11; G4: #8; G5: #1; G6: #7; G7: #6 and #9; and G8: #2, #3, #4 and #5 (Table S1) based on 
successive exclusion from pooled data. 


Table 2 Analysis of variance for the allometric model test 


SV df SS MS F P 
Parameter 24 21,084.4846 - 
Reduction (B's) 2 21,083.7925 - 
Reduction (Ho) 22 0.6921 0.03146 43.18 <0.00001* 
Residual 5925 4.3166 0.00073 
Total 5949 21,088.8011 


Note: Both dependent and independent variables were log-transformed before analysis. Data were derived from the calibration data set 
(n=5950). SV, source of variation; df, degree of freedom; SS, sum of squares; MS, mean squares; P, probability. Reduction (fB's), test if 
all individual models could be explained by a single equation; Reduction (Ho), difference in the number of parameters. * means that 


there is a high probability that, at least, one population differs to the others. "-" means no value. 


3.2 Leaf morphological trait 


Leaf morphology traits revealed significantly differences between the semi-arid and arid 
populations and those from other provenances (Table 3). Our results indicated that the population 
from Paulista-PE (#3) displayed higher leaf width, leaf length and leaf area than the other 
populations (Table 3). Morphological traits displayed highly significant differences among four 
populations of C. procera from Arid region, Atlantic Rainforest, Caribbean and Mountain 
Rainforest (Table 3; Fig. 2). The L/W ratio indicated that leaves from Oaxaca-MX (#12) were 
broader (1.28 (+0.14)) than leaves from Afranio-PE (#10; 2.19 (+0.44)) (Table 3; Fig. 2). The 
SLA value in leaves from Afranio-PE plants were 47% higher than that of leaves from 
Oaxaca-MX. Hence, the relationships between L/W ratio and SLA and between L/W ratio and 
mean annual precipitation showed a negative correlation (Fig. 3). In addition, a positive 
correlation between SLA and mean annual precipitation was found (Fig. 3). 


3.3 Leaf anatomical trait 


C. procera leaves had a uniseriate epidermis. The epidermis was formed by cells with curved 
anticlinal walls, and leaves were covered by a thick cuticle along their entire length (Fig. 4). The 
leaves of all morphotypes were characterized as amphistomatic, with stomata arranged in a plane 
above the epidermis (data not shown). Except the population Afranio-PE (#10), C. procera 
presented a dorsiventral mesophyll parenchyma, consisting of three to four layers of palisade 
parenchyma with slightly elongated and juxtaposed cells that were thin and compactly arranged. 


Table 3 Leaf width (W), leaf length (L), leaf area (LA), L/W ratio and specific leaf area (SLA) of Calotropis 
procera from 12 populations 


W (cm) L(cm) LA (cm?) L/W ratio SLA 
(m?/kg)” 


Pop. No. 
Mean* Min. Max. Mean Min. Max. Mean Min. Max. Mean Min. Max. 


#01 500 11.840.128 1.85 13.54 21.040.14° 2.51 22.17 197.742.97' 4.78 242.87 1.8340.17 1.06 1.83  12.68+0.32%° 
#02 500 15.540.144 0.64 17.29 24.4+1.19 1.13 26.15 309.546.08° 0.55 377.59 1.6940.16 1.35 3.05 9.89+0.9 1° 
#03 499 18.340.21* 0.33 19.97 27.2+0.27* 0.70 32.20 418.6+7.92* 0.20 526.86 1.7340.23 1.15 2.51 10.67+0.56% 
#04 500 16.340.22° 0.49 18.46 25.2+0.24> 0.98 27.90 346.346.98> 0.35 421.16 1.7440.23 1.32 2.68  11.66+0.14¢ 
#05 500 14.5+0.20° 0.29 16.93 23.2+0.36° 0.69 27.68 280.347.664 0.18 382.94 1.8540.26 1.19 2.78 9.98+0.45°f 
#06 499 14.6+0.20° 0.48 16.69 23.5+0.18° 1.04 24.77 275.4+5.654 0.37 331.28 1.76+0.21 1.36 2.62 10.26+1.17%f 
#07 500 14.640.19 1.01 17.12 22.8+0.20° 2.17 24.75 275.946.034 1.53 363.71 1.70+0.20 1.24 2.63 8.89+0.04°f 
#08 491 14.3+0.16° 0.38 16.53 22.340.164 0.95 23.70 252.744.20° 0.32 298.15 1.81+0.28 0.52 2.52 9.58+0.75% 
#09 500 14.940.18° 1.00 16.95 24.3+40.31° 2.23 27.12 295.447.654 1.77 360.02 1.6940.16 1.36 2.26 9.3340.48% 
#10 466 13.240.16' 0.28 14.77  25.040.30° 0.72 27.59 249.545.22° 0.24 294.80 2.1940.44 1.70 6.23 8.5741.36! 


#11 499 = 17.1+0.21 0.70 19.40 24.8+0.32° 0.91 27.50 336.948.79° 0.56 426.75 1.4540.15 1.01 2.38  13.81+0.56° 


#12 493 18.8+0.16° 0.64 20.70 25.1+40.26" 0.90 28.50 351.547.46> 0.51 423.29 1.2840.14 0.21 1.62 16.11+0.53* 


Note: Pop., population code. No., number. * means of 5% bigger expanded leaves (n=25). **, n=50. Mean+SE. Different lowercase letters within a 
column indicate significant differences of means of W, L, LA, L/W ratio and SLA among different populations (P<0.05). 
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Fig. 2 Leaf morphological traits of Calotropis procera populations from the Afrânio-PE (Brazil; a), Paulista-PE 
(Brazil; b), Oaxaca-EUM (Mexico; c) and Montería-CO (Colômbia; d) 
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Fig. 3 Correlation between specific leaf area (SLA) and L/W ratio (L, length; W, width) (a), correlation between 
SLA and mean annual precipitation (b), and correlation between L/W ratio and mean annual precipitation (c) 
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Fig. 4 Leaf anatomical traits of cross-sections of Calotropis procera populations from the Afrânio-PE (a), 
Paulista-PE (b), Oaxaca-MX (c) and Montería-CO (d). ADAE, adaxial epidermis surface; PP, palisade 
parenchyma; SP, spongy parenchyma; ABAE, abaxial epidermis surface. Air spaces in spongy are shown. 
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However, the leaves in Afrânio-PE displayed an isobilateral mesophyll (Fig. 4). In all populations, 
we verified that spongy parenchyma was consisted of loosely arranged cells of various shapes, 


Fig. 5 Leaf anatomical traits of paradermal sections of Calotropis procera populations from the Afrânio-PE (a), 
Paulista-PE (b), Oaxaca-MX (c) and Montería-CO (d). Air spaces in palisade are shown. 


An in-depth analysis was conducted only among 4 populations (#3, #10, #11 and #12). C. 
procera leaves from Afrânio-PE presented thicknesses of a leaf, adaxial and abaxial epidermis, 
and palisade and spongy parenchyma 23%, 26%, 13% and 26% thicker, respectively, than leaves 
from Montería (Table 4). Despite these differences, leaves from Afrânio-PE displayed 3.7 (+0.1) 
adaxial palisade layers, whereas leaves from Montería-CO displayed 3.3 (+0.1) adaxial palisade 
layers. The number of palisade layers was, thus, 13% higher in leaves from Afrânio-PE than in 
leaves from Montería. The palisade and spongy parenchyma thicknesses of leaves from 
Afrânio-PE were 214% and 40% thicker, respectively, than leaves from Montería-CO (Table 4). 


Table 4 Leaf anatomical traits of Calotropis procera from 4 populations 


Population 
Parameter 
Afrânio-PE (#10) Paulista-PE (#3) Oaxaca-MX (#12) Montería-CO (#11) 
Total leaf thickness (um) 792.26+23.46* 769.10+27.70* 689.74+22.05° 643.10412.91° 
Adaxial epidermis thickness (um) 42.10+0.49* 38.08+0.49° 34.47+0.63° 33.46+0.70° 
Abaxial epidermis thickness (um) 34.79+0.47* 34.08+0.87* 30.68+0.58° 30.90+0.61° 
Palisade parenchyma thickness (um) 216.64+4.85° 211.16+5.29° 187.23+3.92° 172.44+5.79° 
Spongy parenchyma thickness (um) 501.82+8.94* 509.58+21.36" 428.99+13.03° 399.73+10.79° 
Adaxial palisade layers 3.66+0.07" 3.62+0.07" 3.42+0.05° 3.32+0.06° 
Air spaces in palisade (%) 12.3940.84° 14.42+1.21° 25.17+1.35° 38.95+1.09* 
Air spaces in spongy (%) 44.58+1.56° 45.35+2.46° 56.8741.69° 62.26+0.64* 


Note: Different lowercase letters within a row indicate significant differences among different populations at P<0.001 level. Mean+SE; 


n=20. 
3.4 Leaf symmetry trait and plasticity index 


The C. procera leaves demonstrated a fluctuating asymmetry (data not shown). The data on the 
plasticity index indicated that, although foliar plasticity index was quite high for some 
populations and lower for others, this index did not significantly differ among the assessed 
populations, i.e., plasticity index did not appear to be provenance-dependent. On the other hand, 
plasticity index indicated that morphological traits (0.77 (+0.09)) showed a higher plasticity than 
anatomical traits (0.34 (+0.04)) (Table S2). 
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3.5 Pearson's correlation coefficient and PCA 


The correlation between mean annual precipitation and leaf morphology was strongly and 
negatively significant based on the L/W ratio (r= —0.41), total leaf thickness (r= —0.31), adaxial 
epidermis thickness (r= —0.35), abaxial epidermis thickness (r= —0.54) and palisade thickness (r= 
—0.32), but positively significant based on the elevation above sea level (r=0.26), SLA (r=0.55), L 
(r=0.28), W (r=0.33), LA (r=0.29) and air spaces in spongy mesophyll (7=0.28; Fig. 6). 
Meanwhile, no correlation between mean annual precipitation and leaf symmetry was observed. 
Across all populations, the effect of elevation above sea level on anatomical characteristics was 
marginal compared to the effect of mean annual precipitation, with no significant correlation 
(P=0.985). 

Significant patterns regarding the morpho-anatomical traits were observed in different 
provenances (i.e., Atlantic Rainforest, Semi-arid region, Arid region, Caribbean and Mountain 
Rainforest). Specifically, for the morphological traits (i.e., L, W, LA and L/W ratio), the first two 
axes of the PCA (cumulative proportion=98.11%) indicated a partial separation of the natural 
morphological variation spectrum in leaves. The leaves collected from the Rainforests and 
Caribbean environments showed greater expression and dominance of these morphological traits, 
as indicated by their relatively larger leaves (Fig. 6). Although overlaps between groups in 
morphological traits were observed, the analysis of variance confirmed the consistent differences 
among populations (P<0.001) except between Atlantic Rainforest and Caribbean that 
systematically presented more similar leaf morphological traits. The anatomical traits also 
indicated a clear separation in different provenances. The two axes explained 66.00% of the data 
variance. It should be noted that the air spaces in the palisade and spongy parenchyma (visible in 
the upper left-hand portion of the PC1 axis) revealed expressive values for the semi-arid more 
humid environments (Rainforests and Caribbean). The other end of the axis was dominated by the 
traits of palisade thickness, total leaf thickness and spongy thickness, representing a specialization 
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Fig. 6 Principal component analysis (PCA) ordination of Calotropis procera from 12 populations ranked on the 
basis of 2 ecological, 6 morphological and 9 anatomical traits. The labels indicate traits with the highest 
eigenvector values on PCA axes 1 and 2; the labels with the highest scores are presented nearest to the axes. The 
stylized figures indicate the extremes of precipitation and leaf morphologies. Leaves from Oaxaca-MX are at the 
higher end of PCA axis 2, and leaves from Afranio-PE are at the lower end of PCA axis 1. L/W ratio, ratio of leaf 
length to leaf width; ADAE, adaxial epidermis thickness; ABAE, abaxial epidermis thickness; PN, number of 
cells layers in palisade mesophyll; AIR, absolute inter-rib value; LA, leaf area; L, leaf length; W, leaf width; ASS, 
air spaces in spongy mesophyll; ASP, air spaces in palisade mesophyll; AMR, absolute midrib value; SLA, 
specific leaf area; P, precipitation; ST, spongy parenchyma; TT, total leaf thickness; PT, palisade parenchyma 
thickness; ELV, elevation above sea level. 


for the semi-arid and arid populations. The anatomical traits showed significant differences 
among different provenances (P<0.001; Fig. 6), except between Arid region and Atlantic 
Rainforest. A coordinated change in the two-dimensional axis indicated that an increase in the air 
spaces of palisade mesophyll was related to a decrease in the thickness of the adaxial epidermis. 
Meanwhile, differences in the climatic variations were significant associated with the 
morpho-anatomical changes (P<0.001). Similarly, the SLA also showed a precipitation-related 
change. 


4 Discussion 


This study demonstrated that C. procera plants collected from wet-to-dry environments have 
developed a wide range of morphological traits, resulting in unique model (for each morphotype 
or small population) for leaf area estimation based on linear leaf dimension. In addition, 
anatomical and physiological adaptations in their leaves are also presented. Thus, the major 
differences found in the L/W ratio and the SLA in C. procera populations may be the most 
important traits for the survival in dry environment. Parkhurst and Loucks (1972) state that 
"natural selection leads to organisms having a combination of form and function optimal (i.e., 
phenotypic plasticity) for growth and reproduction in the environments in which they live". Thus, 
through natural selection, the form and function of the members of a population may be expected 
to approach or reach an optimum performance. In this case, C. procera presented distinct models. 
This behaviour is divergent from Coffea spp. (Antunes et al., 2008) or Jatropha curcas trees 
growing in wet or dry environments (Pompelli et al., 2012). One can ask: why leaves are narrower 
in dry environments and broader in wet ones? Maximum leaf size undergoes a considerable 
acclimation to water stress, whereas leaf shape is largely controlled by heredity (Arterburn et al., 
2010). Plant ecologists have emphasized the extensive relationships between leaf traits of native 
plants and the climate (Wright et al., 2004; Dainese et al., 2015; Diaz et al., 2016). For example, a 
common strategy in plants inhabiting arid and semi-arid regions is the development of thicker 
leaves with low SLAs (Bacelar et al., 2007; Evert, 2013). However, the coordination between 
foliar traits in the multivariate environments seems to be a more complex and, at least partially, 
climate-dependent mechanism (Wright et al., 2004). 

Rivas et al. (2017) found that temperature of leaves of C. procera growing in wet environment 
are lower than that of air (-1.33°C), whereas temperature of leaves growing under drought are 
12.61°C higher than that of air. Under drought, C. procera reduces stomatal conductance and 
increases leaf temperature (Antunes et al., 2017; Rivas et al., 2017). The high plasticity of leaf 
would contribute to radiation interception and radiative and sensible heat loss, increasing the 
intrinsic water use efficiency of C. procera that can partially mitigate the extreme water vapor 
gradient between mesophyll and the surrounding environments. Similar results are described in 
Quercus coccifera under two Mediterranean environments (Peguero-Pina et al., 2016b). 
Multi-layered parenchyma have been identified in a variety of plant families (Fahn, 1990; Evert, 
2013), however, variation in parenchyma thickness and layers in response to a natural gradient in 
precipitation is firstly reported. Based on this, we argue that more palisade layers or an 
isobilateral mesophyll may protect the photosynthetically tissue (Chazdon and Kaufmann, 1993). 
Michaletz et al. (2015) found that lowering transpiration by altering leaf morphology would 
reduce the loss of hydraulic conductance by xylem embolism under drought. Larger leaves 
present smaller convection coefficients and higher values of resistance to heat transfer through 
leaf boundary layers than smaller leaves, since leaf size may change for optimizing leaf 
temperature (Leigh et al., 2016). Narrower, smaller and thicker leaves in drier environments 
exhibit greater loss of radiative and sensible heat, with lower latent (transpiration) heat losses 
(Peguero-Pina et al., 2014). Thicker leaves retain more water and become warmer, mainly due to 
their high absorption of near-infrared radiation by water, which dissipate energy absorbed via a 
non-transpiring system. They open their stomata less frequently and close them early in the day 
(Rivas et al., 2017) for saving water and displaying a drought postponement morphotype 
(Antunes et al., 2017). In this way, low SLA is important trait for the adjustment of leaf energy 
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budget, significantly increasing radiative heat loss in dry environments. Moreover, sclerophylly 
(low SLA) is very common in species that grow in semi-arid and arid regions, and this trait may 
function to save-water (Monk, 1966). C. procera present a morphotype with thinner leaves in wet 
environments, while thicker leaves in dry environments. Thinner leaves dissipate incident energy 
at the expense of a higher rate of transpiration (latent heat loss). Thinner leaves with larger areas 
may lead to a physiological consequence of a high transpiration rate of entire plant. However, it 
should be considered as a neutral trait for plant survival because water is available in plentiful 
amounts in their original environment. When we associate this trait with higher stomatal 
conductance, increased CO> influx, larger leaf area and greater net photosynthesis, thinner leaves 
of C. procera is a positive trait. The thinner leaves of C. procera allow it to survive in dry and wet 
environments. These findings are consistent with Lindorf et al. (1991) and Witkowski and Lamont 
(1991), who previously reported that C. procera leaves are more succulent in arid environments, 
with a thicker cuticle and small leaf area. 

Thinner leaves are often short-lived (e.g., those of plants from Monterfa-CO (#11) or 
Oaxaca-MX (#12)), which may be associated with lower leaf area construction costs (Villar et al., 
2013). Lower area-based leaf construction costs were verified for invasive Lythrum salicaria 
plants compared with five co-occurring herbaceous plants and shrubs (Nagel and Griffin, 2001). 
The authors suggested that the production of a larger photosynthetic leaf surface area at a lower 
energy expense, resulting from either lower requirements or higher energy efficiency, would help 
the plant to outcompete neighbouring native species in low stress, low disturbance and 
high-competitive wet environments. Research found that increased thickness of the mesophyll or 
a more compact mesophyll would improve CO2 diffusion through the intercellular spaces from 
the sub-stomatal cavity to the outer surfaces of the mesophyll cells (Tezara et al., 2011). In the 
present study, we imply that plants in more xerophytic environments might have lower stomatal 
conductance to mitigate possible water losses to the atmosphere, without large losses of 
photosynthetic capacity, thus leading to an increase of WUE;, as previously reported by Rivas et al. 
(2017). 

Leaf area changes systematically with leaf dimensions. The morphotype from Afranio-PE 
appeared better adapted to drought with a smaller leaf. The result indicates that selecting thinner 
and larger leaves, as well as lower L/W ratio is associated with maximum light capture at the 
expense of abundant water use (White and Montes-R, 2005). A thicker palisade parenchyma, 
higher number of adaxial palisade layers and less air space may exert important influence on the 
photosynthetic capacity of the leaf. Thicker leaves with lower SLA present a greater amount of 
living photosynthetic cells per unit of leaf area (Villar et al., 2013). In arid environments, leaves 
are often submitted to very high direct irradiance and great fluency of diffuse light. Thus, thicker 
leaves absorb more light than thinner ones (Falcioni et al., 2017). Upper palisade parenchyma 
cells are very efficient in absorbing blue and red lights (Terashima et al., 2009; Falcioni et al., 
2017; Smith et al., 2017). However, sunlight is rich in yellow and green lights (Terashima et al., 
2009) that are weakly absorbed and have been proposed to penetrate deeper in leaves, allowing 
for light reactions in the middle spongy parenchyma. Isobilateral symmetry of palisade 
parenchyma found in the Afranio-PE may be another adaptive trait to absorb soil-reflective and 
diffuse lights, efficiently driving photochemical reactions of photosynthesis on both sides of the 
leaf lamina, and reducing the pressure for transpiration caused by larger leaf areas. 


5 Conclusions 


C. procera leaves demonstrate a considerable capacity to adjust their morphological, anatomical, 
and physiological traits to different environments. C. procera presents a great plasticity, with the 
values ranging from 0.115 and 0.999. Moreover, morphological traits of C. procera are more 
plastic than anatomical traits, indicating that leaf morphological trait of C. procera is profoundly 
affected by environmental conditions. For example, increase in leaf thickness, lower L/W ratio 
and isobilateral mesophyll in Afranio-PE morphotype strongly contributed to the adaptation to dry 
environments. In this study, at least three quarters of morpho-anatomic traits were found in 


different environments, which was consistent with the change and efficiency of C. procera 
eco-physiological mechanisms. We demonstrated that morphological traits are good predictors for 
the form and function of plants at macroscale. 
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Fig. S1 Deviation of the estimated leaf area from the observed leaf area of Calotropis procera from different 
groups of population. The group information of population is listed in Table S1. Leaf area of C. procera was 
estimated using four models in which Bo and Bı are coefficients. L, length; W, width; e; is the random error that 
represents the discrepancy in the approximation. Vertical bars denote mean value and spreads denote 95% 
confidence intervals of the difference. Asterisks denote significant difference (P<0.05, Student's ¢ test) from the 
observed values. See Figure 1 and Table 1 for details of each population. 


Table S1 Group, coefficient (Bo), standard error (SE), coefficient of determination (R?), degree of freedom (df), 
sum of square (SS), calculated F, P and model of estimated leaf area as a function of linear dimensions of leaves 
(L, length; W, width) of Calotropis procera 


Group Coefficient SE R? df SS F P Estimate of leaf area (Y) 
(---)/All Bo=0.71924 0.00808 0.9894 5945 521,536 557,318  <0.0001 Y=0.71924x(Lx W)!-01867 
(G1)/#12 Bo=1.22999 0.11754 0.9040 491 325,645 4632 <0.0001_ = Y=1.22999x(Lx W)991999 
(G2)/#10 Bo=0.75486 0.01925 0.9959 464 9151 114,472 <0.0001 Y=0.75486x(Lx W)!.00023 
(G3)/#11 Bo=0.61262 0.01294 0.9980 497 10,010 256,617 <0.0001 Y=0.61262x(LxW)!04114 
(G4)/#8 Bo=0.66429 0.01072 0.9990 489 2763 488,336 <0.0001 == Y=0.66429x(Lx W)!.03074 
(G5)/#1 Bo=0.64241 0.00981 0.9976 498 3134 206,263 <0.0001 Y=0.64241x(LxW)!.03857 
(G6)/#7 Bo=0.69267 0.01140 0.9974 498 6661 192,771 <0.0001 Y=0.69267x(LxW)!.931%6 

(G7)/#6, #9 Bo=0.68226 0.00866 0.9971 997 16,503 343,227 <0.0001 Y=0.68226x(LxW)!0?95° 
(G8)/#2, #3, #4, #5 Bo=0.63834 0.00668 0.9971 1997 57,371 686,867 <0.0001 Y=0.63834x(LxW)!-0426° 


Note: #01, Aimorés-MG; #02, Jaboatão Guararapes-PE; #03, Paulista-PE; #04, Pombos-PE; #05, Santa Terezinha-PB; #06, 
Alagoinha-PE; #07, Jaramataia-AL; #08, Serra Talhada-PE; #09, Parnamirim-PE; #10, Afrânio-PE; #11, Montería-CO; #12, 
Oaxaca-MX. 


Table S2 Plasticity index of leaf morphological and anatomical traits of Calotropis procera from 4 
populations 


Afranio-PE Paulista-PE Oaxaca-MX Monteria-CO 


Morphological trait (#10) (#3) (#12) (#11) Mean+SD 

L/W 0.79 0.54 0.87 0.58 0.68+0.06 

Specific leaf area 0.31 0.12 0.22 0.24 0.22+0.04 

Leaf length (L) 0.97 0.98 0.97 0.97 0.97+0.02 

Leaf width (W) 0.98 0.98 0.97 0.96 0.97+0.01 

Total leaf area 0.99 0.99 0.99 0.99 0.99+0.01 

Absolute value between midribs 0.76 0.86 0.74 0.84 0.80+0.02 

Relative value between midribs 0.79 0.86 0.71 0.84 0.80+0.03 

Absolute value between interribs 0.80 0.86 0.87 0.92 0.86+0.02 

Relative value between interribs 0.80 0.88 0.90 0.93 0.88+0.02 

Mean+SD 0.79+0.10 0.79+0.05 0.81+0.08 0.81+0.08 0.77+0.09 
Anatomical trait 

Total leaf thickness 0.32 0.38 0.45 0.28 0.36+0.03 

Abaxial epidermis thickness 0.17 0.35 0.27 0.31 0.28+0.03 

Adaxial epidermis thickness 0.16 0.18 0.30 0.29 0.23+40.03 

Palisade mesophyll thickness 0.30 0.30 0.28 0.42 0.32+0.03 

Spongy mesophyll thickness 0.21 0.39 0.43 0.34 0.34+0.04 

Palisade number 0.24 0.25 0.19 0.24 0.23+0.01 

Air spaces in palisade 0.59 0.71 0.62 0.38 0.57+0.06 

Air spaces in spongy 0.40 0.60 0.27 0.15 0.36+0.08 

Mean+SD 0.34+0.06" 0.30+0.05* 0.30+0.03* 0.35+0.05" 0.34+0.04 


Note: * denotes statistically significant differences among anatomical traits in different populations at P<0.01 level. Mean values for the 
leaves collected from different populations among morphological traits did not differ significantly at P<0.05 level. 


